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Abstract 

For a Friedman-Robertson- Walker space-time in which the only 
contribution to the stress-energy tensor comes from the renormalised 
zero-point energy (i.e. the Casimir energy) of the fundamental fields 
the evolution of the universe (the scale factor) depends upon whether 
the universe is open, flat or closed and upon which fundamental fields 
inhabit the space-time. 

We calculate this "Casimir effect" using the heat kernel method, and 
the calculation is thus non-perturbative. We treat fields of spin 0, ^, 1 
coupled to the gravitational background only. The heat kernels and/or 
zeta-functions for the various spins are related to that of a non- 
minimally coupled scalar field which is again related to that of the 
minimally coupled one. A WKB approximation is used in obtaining 
the radial part of that heat kernel. 

The simulations of the resulting equations of motion seem to exclude 
the possibility of a closed universe, K = +1, as these turn out to have 
an overwhelming tendency towards a fast collapse - the details such 
as the rate of this collapse depends on the structure of the underly- 
ing quantum degrees of freedom; a non-minimal coupling to curvature 
accelerates the process. Only K = —1 and K = will in general 
lead to macroscopic universes, and of these K = —1 seems to be more 
favourable. 
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The possibility of the scale factor being a concave rather than a con- 
vex function, as is the case for a K = —1 FRW-spacetime inhabited 
by a conformally coupled scalar field, potentially indicates that the 
problem of the large Hubble constant is non-existent as the age of the 
universe need not be less than or equal to the Hubble time. 
Note should be given to the fact, however, that we are not able to pur- 
sue the numerical study to really large times neither to do simulations 
for a full standard model. 

1 Introduction 

In a recent paper we realized that the Casimir effect may drive infla- 
tion or otherwise strongly influence the evolution of space-time because the 
vacuum fluctuations couple to the (background) gravitational field. As one 
cannot escape the Casimir effect it may thus be that this property of the 
zero-point energy can be used for determining whether the universe is open, 
flat or closed. Or working the other way around, one may be able to use 
cosmological observations to tell us about the particle content in the very 
early universe. 

In this paper the quantum flelds reside in a Friedman-Robertson- Walker 
space-time. In section 2 we show how one in general calculates the Casimir 
effect of a free scalar fleld using the heat kernel method and continue to deter- 
mine it explicitly in section 3 in the case of the Friedman-Robertson- Walker 
metric and it is noted that using this approach, the zero-point energy is al- 
most automatically renormalised. The renormalised Casimir energy is then 
inserted in section 3.1 into the right hand side of the Einstein equations (it 
gives an effective action for the scale factor a{t), the only degree of freedom) in 
order to obtain the equations of motion and the results of the corresponding 
simulations of the evolution of the space-time, perhaps relevant if the Higgs 
fleld is believed to dominate the early universe, are presented in section 3.2. 
Having established a procedure for determining the Casimir effect on the 
evolution of space-time we continue to consider the case of a non-minimally 
coupled scalar fleld (in section 4) and show how the Casimir energy is related 
to that of the minimally coupled scalar fleld. Then we proceed in sections 
4.1 and 4.2 to present the equations of motion and the corresponding simu- 
lations for a conformally coupled scalar fleld. In section 5 the approach used 
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in section 4 is generalised to the case of vector bosons, with the problematic 
self-interaction being treated using a mean field approximation procedure - 
due to the complexity of the solution we have not been able to carry out the 
simulations, however. Finally in section 6, after some general considerations 
about how to link a fermionic operator to a scalar or bosonic one, we take 
on free (spin 1/2) fermions. For completeness, we have included a discussion 
of fermions interacting with gauge fields in an appendix (based on the mean 
field approach developed for pure Yang-Mills theory). Einstein equations and 
simulations for free fermions are presented in sections 6.1 and 6.2 Section 7 
gives a few concluding remarks. 

2 General Considerations concerning the cal- 
culation of the Casimir Effect for a Free 
Scalar Field 

We need to determine the right hand side of the Einstein equations, the 
energy-momentum tensor. Now, the quantum expectation value, {T^^), (the 
zero-point energy) is given in terms of the effective action F as [0, § 

(V) = 2(-9)-'^^ (1) 
which is just the quantum analogue of the classical relationship 

jiclass _ _„-l/2 



Another way of looking at this (the way chosen in |l[]), is to use an analogy 
with classical statistical mechanics and introduce Helmholtz' free energy, F, 
which is related to the partition function by 

F = -ilnZ (2) 

where (3 is the inverse temperature and where the partition function is given 
by the functional integral 

Z = J e-^Vcj) = e-^^ff (3) 
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with Feff being the effective action, i.e. 



T = PF = - \nZ (4) 

Actually, this is an effective potential for the gravitational degrees of freedom, 
the full effective action for gravity being given by 

T,^n[g] = SEn[g] + r[g] (5) 

but we will always take the Einstein-Hilbert action S'eh to be understood, 
and simply refer to F as the effective action. This slight abuse of notation 
should not cause any misunderstandings. 

IFtoiu. Helmholtz' free energy one can then, as in [^, proceed to find the 
internal energy and the pressures. It is because of this analogy with the 
free energy we refer to the entire scenario as Casimir driven evolution, the 
renormalised zero point energy being analogous to the well known Casimir 
effect in flat space. 

Since the effective action is given by integrating out the matter degrees of 
freedom in the partition function, it gives the effective equation of motion 
for the remaining degree of freedom, namely the gravitational one, which for 
a Friedman-Robertson- Walker background reduces to an effective equation 
of motion for the scale factor a{t). 
Now, since we're dealing with a free particle, 

S=^l d^x^ig^^d.^d,^ - (6) 

albeit in a special geometry, this integral is merely a Gaussian. Thus we can 
perform it using standard techniques, arriving at the result 

Z = (detA)-^/2 
where the differential operator A is given by 

^.g^n-M^ (8) 
with □, the d'Alembertian of the curved space-time, given by ^ 

□ = -^d, (g^''J\^\^^ (9) 
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The effective action is tlius seen to be given by a functional determinant; 

r = -ln(det(n-M2))"^^^ (10) 

Tlie determinant of an operator, A, is, by definition, tlie product of its eigen- 
values 

det{A)=Y[X (11) 

A 

The zeta-function is given by 0, Q 

a(5)=E^"' = Tr A- (12) 

A 

and it is related to det{A) by the following equation 

dCA 



ds 



= ^-InA- A-'|,=o = - J^lnA = -InJlA = -lndet(A) (13) 

OA A A 



Now (as shown below) the zeta-function, (a, can be determined from the heat 
kernel Ga{x, x', a) which is a function satisfying the heat kernel equation^ 

d 

A^Ga{x, X, a) = -—Ga{,x, x\ a) (14) 

subject to the boundary condition Ga{x^x\Q) = 6{x,x'). 

To see how Ga{x,x' ,a) is related to (a note that if A(j)x = X4>x then 

GA{x,x',a) = J2x4>x{x)(f)*x{x')e~^'" satisfies the heat kernel equation. This 

in turn makes it possible to show the following relationship between Ga and 

Ca 

roo r roc r 

I daa'-^ I y/^d^xGAix,x,a) = I daa^^W v^d^x ^ |0A(a;)|^e-^" 

^ A 

roo r 

= EX daa^-'e-^'^ J ^d'x\M^)\' 



^If A is the Laplace operator and the space is then becomes (proportional to) 
the Riemann zeta function, and if the space is and A — the heat kernel equation 
becomes the usual heat equation. 
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where, with suitable normahsation, the last integral can be put equal to unity 
which gives us 



poo r roc 

J daa^~^ J ^d'^xGA{x,x,a) = daa 

roo 

= y\^' di\a){\a 
^ Jo 



A 



n-sKAis) (15) 



so that 



Ca{s) = J daa"" ^ J GA{x,x,a)y/gd^x 



and so, finally, for the minimally coupled scalar case 

If d r°° (7'^~^ 

^=2] ^^'^ Ts Jo d^Y{rf--m<^^^^^) (16) 

Note that the integral over x is taken along the diagonal x = x'. Also note 
that the cr-integral is a standard one if one uses the spectral representation 
of the zeta-function, thus easy to calculate though it generally needs further 
regularisation (in the case of a minimally coupled field we use the spectral 
representation of the heat kernel (in a WKB approximation) while in other 
cases we use a hybrid expression, partially referring to eigenfunctions, partly 
a Taylor series generalising the usual Schwinger-DeWitt expansion |T^), in- 



cluding the exponential of the scalar curvature as shown by Jack and Parker, 



[0, to be needed. 

Thus to calculate the zero-point energy of a scalar field one determines the 
scalar field operator (the d'Alembertian) in the relevant space-time, solve 
the corresponding heat kernel equation and subsequently calculates the zeta- 
function from which the generating functional and thus all relevant quantities 
(including the zero-point energy) can easily be calculated. 
The zeta-function of spin 1 gauge bosons can be determined by a generalisa- 
tion of the procedure for the non-minimal scalar field while the zeta-function 
associated with fermions is directly related to that of the a non-minimally 
coupled scalar field. This in turn is related to the zeta-function of a minimally 
coupled scalar field, substantially reducing the amount of work involved in 
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determining the zero-point energy of these fields, as will be shown in the 
following sections. 



3 Solving the Heat Kernel Equation for the 
Minimally Coupled Scalar Field 

Having in principle established a procedure for determining the zero-point 
energy for a minimally coupled scalar field the next step is to find the heat 
kernel and subsequently the Casimir energy of the field explicitly. The heat 
kernels for higher spin, as well as for non-minimally coupled scalar fields, are 
going to be related to the result found in this section. Thus this calculation 
is a central part of the paper. 

If one writes the Friedman-Robertson- Walker line-element with the scale fac- 
tor a — a{t) as 

ds'' = dt'' - a^{t){dx^ + f{x){de^ + sm\e)d(l)^)) (17) 

where 

sinx for = (closed universe) 
f{x) = {x ior K — (fiat universe) (18) 
sinhx ior K — —1 (open universe) 

then the d'Alembertian becomes 

□ = -J=a^ (<,<"-/^s„^ (19) 

where L is the usual angular momentum operator. 

Even though the equation separates, the minimally coupled scalar heat kernel 
equation 



nGu{x,x\a) = {{d^ + 3^dt)-^Jdi-2d^{^d,)]-^^]Ga{x,x',a) 



~Gu{x,x',a) (21) 
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is rather hard to solve. Fortunately, the eigenfunctions, ipx{x), of the spatial 
part of the d'Alembertian are already known to be Q 



(27r) 2 6^^'^ where k = {ki, k2, k^) for = 

^uix)Yf{e, (j)) where k = (fc, J, M) for K = ±1 



where 

— oo < fcj < oo ; where k = |k| for = 
M = - J, -J+ 1,..., J; 



J = 0, 1, A; - 1 ; k= l,2,...foT K = 1 
J = 0, 1, ... ; < A; < oo for i^' = 



and 

i+J 

d cosh X 



Ul~j\x) = [l7rk\e + l)...{e + {2 + jY)^ %inhx(^^^) ' cos(A:x) 



( d \^^^ 

xisinxH cos(fcx) (24) 

\ " cos X J 

as quoted by so proceed by assuming the heat kernel of the d'Alembertian 
to be of the form 

Gn=^7AA(x)^:(x')rA(t,t',a) (25) 

A 

which is a solution of the heat kernel equation provided that the set of func- 
tions Tx{t,t',a) solves 

(^2 + S-dt)Tx{t, t', a) = -d^t, t', a) (26) 
a 

Writing 

Ut,t',a) = a-'/\t)Mt)a-'/'{t')Mt')e-^- (27) 



gives the following equation for the function /, 



A 



27 /a\2 3a 



8 



with 

X = k'^-K (29) 

being the eigenvalue of the spatial part of the d'Alembertian. 
We have not been able to solve this equation explicitly (due to the fact that 
we have no a priori knowledge about hx{t)). Note however the similarity with 
the stationary Schrodinger equation; ^[E — V{x)]ip{x) = /\ip{x). Thus if 

^ is negligible it is sound to use the WKB-approximation0 giving 

/, = /i7(t)e/o^*'v^(*') (30) 

One should note that this is the only approximation introduced so far. 
Having obtained a set of functions on which to expand the heat kernel; 

Gu{x,x',a) = 5]0A(a;)0I(a;')e-^-'^ 

A 

= 5:^A(x)V^I(xOrA(t,t',a) 

A 

= ^^,(x)^Kx')«-'/'W/A(t)a-'/^(t')/A*(i')e-"'^ 

A 

= E^A(x)^I(x')a-^/2(t)/i7(t)e/o'^*"^(*") X 

A 

a~''/\t')hl^{t')Jo '''"'^^^'"'^e-^'' (31) 

the C-function reads 

/ h~'/\t)e~'foV^)d^'e'''^a'-'dtda (32) 



where we have used ^/g = a^f'^ sin 9 = a?\J~g^ and that the eigenfunctions 
of the spatial part of the d'Alembertian are orthonormal with respect to the 

measure \J~g^cPx = p{x)dxd^- 

^This approximation assumes that h and thus the scale factor a is relatively stable and, 
in turn, tends to stabilise a. Thus, when simulating the behaviour of the scale factor later 
on, the evolution of a will probably be understated. 



9 



3.1 The Einstein Equations for the Friedman-Robertson- 
Walker space-time with a minimally coupled scalar 
field 

With the metric given by equation (^) the non-zero components of the 
Einstein tensor are (/' = d^f) 

_ 2//" + r-3aV'-l 
tjoo — 



G 



/'2 + (-2aa - a2)/2 - 1 
/" + (-2aa - a^)f 



22 



_ r + {-2aa - g^)/ 
G33 - (33) 

Inserting / from equation (19) it is exphcitly seen that, incidentally, the 
(two) last components can be rewritten as 

= G,. = 2^i±4!±iL (34) 

leading to the following (fourth component of the) equation of motion for the 
scale factor a{t) (using the first or second component would have complicated 
things unnecessarily as t would no longer be the only variable) 

2ad + a^ + K 

= (^33) (35) 

One should note that while this may sound a bit arbitrary (involving a choice 
of using only G33), it is actually the uniquely defined equation of motion 
which follows from the effective action. The full Einstein tensor is only listed 
in order to give an interpretation of the resulting equation, namely as the 
pressure coupled to G33. 

We proceed to calculate the resulting effective action exphcitly using equa- 
tions (11,13,11): 
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d 

ds 



s=0 



J2^is) Jhxitr'/'e-^Io^'^''dt 



(37) 



with 



hx = 3da 



-2 



3.. 

—aa 
2 



-1 



+ A + a^^A 



(38) 



coming from the WKB-approximation. The right hand side of the equation 
of motion (^) for the scale factor a{t), obtained by varying F with respect 
to a{t), thus reads 



InA /i 



-3/2 



'-aa 
2 



-aa 
4 



Aa 



+ 



^1/2 



-Qda ^ + |aa ^ — 2Aa ^ 



3da ^ — |da^i + A + a^^A 



dt'le-^/o^-^*' (39) 



We note the appearance of an integral over time: Formally the system de- 
pends on its entire past. Even though this dependence makes simulations 
somewhat tedious, we have still carried out quite a few, and these will be 
discussed in the next subsection. 

These equations differ from the ones found by Suen and Anderson, [Q, due 
to the use of the WKB-approximation. Note by the way that it is the expo- 
nential coming from this approximation which ensures rapid convergence of 
the summation over A, thereby making the numerical solution more feasible. 



3.2 Simulating the Equations of Motion 

In order to simulate the evolution of a Friedman- Robertson- Walker spacetime 
according to the above equation of motion for a{t), we first have to discretise 
the time, i.e., to introduce a time-step St. One immediate problem is the 
appearance on the right hand side (the free energy contribution) of a{t) and 
its derivatives raised to various powers and exponentiated. Thus we cannot 
just insert d(t) = {a{t + 6t) -a{t))/6t, d{t) = {a{t + 5t)-2a{t) + a{t-5t))/5t^ 
into the equation of motion and then find a recursion relation telling us how 
to compute a{t + 5t) from the knowledge of a(t), a{t — St), as we are not able 
to isolate a{t + 6t). Instead we will consider the right hand side as a source 
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term to be evaluated att — 6t whereas the left hand side, the Einstein tensor, 
is to be evaluated at time t. One can either justify this by sheer necessity or 
argue that it is plausible that the disturbance caused by the vacuum fluctu- 
ations do not give rise to an instantaneous change of geometry. In the limit 
St ^ this distinction, of course, disappears. With this, we can now isolate 
a{t + 6t) as it only appears in the discretised version of the Einstein tensor. 
As mentioned in the introduction this work was inspired by the simulations 
of the (space-time called a) hyper-spatial tube. Those simulations did show 
a number of qualitatively different types of behaviour, the behaviour being 
determined by the initial conditions. Thus we have varied the initial condi- 
tions which for a{t = 0), in units of the Planck length, were 3,5 or 10 and 
for a{t = 0), in imits of the velocity of light, were 0.0, ±0.1, ±0.5 ± 1.0 and 
1.50. 

Figures la-c shows, for the various initial conditions, simulations for 
K = ±1,0 and —1 respectively. As can be seen, the behaviour of the scale 
factor is rather universal in the sense that for given K the initial conditions 
determine whether the universe will collapse or expand but once it expands 
there is little (qualitative) dependence on the initial conditions of the evolu- 
tion. 

We notice that for = ±1 there is a marked tendency to collapse, whereas 
for the two other cases = — 1, we have a tendency to expansion. Due to 
the amount of computer time and storage required, we have unfortunately 
not been able to pursue these evolutions very far, and so extrapolations on 
the evolution for large t becomes somewhat speculative. However, to see if 
the behaviour was qualitatively stable at somewhat higher times we prepared 
a plot of the evolution for i^' = —1 up to times lO^Tpianck, this did not give 
rise to any new behaviour and is therefore not included here. Also note that 
the expansion becomes polynomial at late timesQ (powers 0.6 — 0.7 for K = 
and ~ 0.8 for K = -1). 

The early paper by Fischetti, Hartle and Hu ||T3| found a solution for K = 



to be of the assymptotic form a{t) ~ which was corrected by Starobinsky 
[IIH] a few years later where he found a(t) ~ 01^^/3(1 + (2/3Mt) sinM(t-ti)), 
which is consistent with our findings. 

■^These are not observable velocities but geometrical ones and thus potentially can be 
larger than c. 

^This is not obvious from the plots, but by ignoring the data from first lOTpjancfc one 
gets data that indeed are on straight lines in a doubly logarithmic plot. 
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It should be noted, however, that both of these papers had to introduce un- 
known constants, ^2,^3 (in the notation of Starobinsky), which they then 
supposed to be both positive (Fischetti et ah) or k2 > 0, < (Starobin- 
sky). These quantities should be determined from the summation over quan- 
tum degrees of freedom. Our approach has no such arbitrary constants since 
we start from first principles, and hence these constants take on definite val- 
ues which apparently are consistent with Starobinsky's results. 
Moreover, both of these early papers take the energy momentum tensor to 
obey p = \Pi which is not necessarily true for quantum corrections, and 
moreover do not violate any energy conditions as the true quantum field the- 
oretical vacuum supposedly does. Since we use the Casimir energy density 
our energy momentum tensor generally violates the classical energy condi- 
tions, as one would expect a quantum vacuum to do. 

The FRW-space-time is somewhat similar to that of the hyperspatial tube [|I| 
(which is a cylinder in 5d, i.e in 4d it has the topology Rx Rx S^). Therefore 
the universal behaviour as well as the slow expansion rate is rather surpris- 
ing: In analogy with the case of the hyperspatial tube (in one of the five 
qualitatively different scenarios), we get expansion by first decreasing the 
size of the universe a bit; at some point, then, a repulsive force sets in, and 
the expansion starts. Looking at the data in a double logarithmic plot shows 
this clearly. But in the case of the hyperspatial tube this repulsive force ul- 
timately gave rise to inflation-like expansions followed by polynomial growth 
(cf figure 5b of Jll), this is not seen here. 

It is of course still possible that such a type of expansion can occur and that 
we simply missed it due to our choice of initial conditions^ (or that some 
other residing field(s) than those considered in this paper could turn on such 
behaviour) . 

4 The Non- Minimally Coupled Case 

The heat kernel equation for a non-minimally coupled scalar field is 

(□ + ^R)Gu+^Rix, x', a) = -d^GuHRi^, a;', a) (40) 

^The calculations done in this paper are 'fully dynamical' while those of 0] are quasi- 
static. However, we would not expect that this is what causes the difference. 
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where R = 6[aa + + K]/a'^ {K = 0, ±1) is the curvature scalar. We know 
how to solve this for ^ = (i.e. minimal coupling), and the solution to the 
generally coupled heat kernel equation is now assumed to be of the form 

Gn+^R{x, x' , a) = Gn{x, x , a)F^{x, x' , a) (41) 

where Gn is the heat kernel for the minimally coupled case (^ = 0). Inserted 
into the heat kernel equation this gives 

(□ + (V In(G'n)) • V + ^R)F^ = ~d„F^ (42) 

where V In Gn ■ V-F^ is short-hand for 

VlnGn ■ VF^ = -^d,{^/g9n{UGnF^) + 2g^^- d^lnGnduFi^ (43) 
Gn^ 

By considering the ^i?-terms as a kind of mass term, one would expect F^ to 
be of the form exp(— ^i?cr + ...). We will write 

F^ = (44) 

and then Taylor expand in powers of a, i.e. 

oo 

^5 = E ^n^" (45) 

n=0 

The boundary condition implies tq = 0, and we can arrive at a recursion 



relation for the remaining coefficients by inserting (p^ and (45) into 
But before we can collect powers of a we must also Taylor expand VlnGn 
as this also depends on a. Thus, we must write 

oo 

V^lnGn = (46) 

n=0 

We then obtain 

-5]n(T"-V„ = 2 5](e;;:Var^)a"+™ + ^(VV„)(Var„)a"+'" + 

n,m n,m 

^(□r>" + ei? (47) 
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i.e. 



'n+l 



n + 1 



. n'=0 



which leads to the following first few coefficients: 



n > 1 
(48) 



To 
T3 







^^g-VaR-leiVR)'- 



and so on. We notice first of all that r„ contains higher and higher deriva- 
tives of i? as n grows. We furthermore notice that, in accordance with our 
expectations, = —^Rcr + Many years ago it was argued by Parker 

and coworkers that such a term should be present in a resummed heat-kernel, 
see e.g. |jT^, but probably due to technical difficulties this suggestion doesn't 
seem to have been taken up by other authors. By exponentiating not just 
the scalar curvature but also the remaining corrections to the heat kernel of 
the d'Alembertian in the way we do it here, we effectively circumvent those 
technical difficulties as will be apparent from what follows below. 
To be able to actually evaluate these coefficients, we need to know Q^. One 
easily sees (by inserting the spectral decomposition of the heat kernel into 
equation (^) and putting cr = 0) that 



^0 



(49) 



where ip\ denotes the eigenfunctions of □. Now, this sum is actually a sum 
over projections when x = x'. As the eigenfunctions are supposedly complete 
this means that the sum equals unity, whereby 



^0=0 



(50) 
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The next coefficient is found by differentiating once with respect to a and 
then putting a to zero. It is 



Again, along the diagonal x = x' the denominator is unity. The numerator is 
actually a spectral decomposition of the d'Alembertian, and hence, inside an 
integral (such as that appearing in the definition of the (^-function), we can 
substitute V^n for (possibly up to boundary terms, which we will simply 
throw away). 

/^From a physical point of view, it is reasonable to only include ffist and 
second derivatives of the scalar curvature R; higher derivatives would be 
higher order quantum effects of the gravitational background, and we are 
only treating a semi-classical approximation to quantum gravity. With higher 
order derivatives of R discarded we thus have 

rs ^ -leiVRf (52) 
r„ ^ n>A (53) 

leading finally to the formula 

Ga+^x, X, a) ^ Gu{x, X, a)e-^^'^+^«°^'^'-^«'(^«)'-' (54) 

A very useful formula indeed, allowing us to express the heat kernel of the 
non-minimally coupled scalar field in terms of that of the much simpler mini- 
mally coupled field. Also, in physical terms, one can think of the ffist term in 
the exponent as due to the background field and the next two terms as due to 
higher order corrections of the gravitational field. Thus, as general relativity 
is only applicable to the 1 loop level, it is probably beyond the theoretical 
framework employed in this paper to expand the exponent to even higher 
order. 

The equation (0) is actually much more general than we need for the 



Friedman-Robertson- Walker geometry, as in this case the curvature depends 
only upon time. In any case, the trick of writing as an exponential and 
the way we were able to eliminate certain derivatives was what allowed us to 
be able to find this result generalising the old result by Jack and Parker [ITl]. 
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4.1 Effective Action and the Einstein Equations in 
Friedman-Robertson- Walker Space-Time for the 
Non-Minimally Coupled Scalar Field 

In order to find the free energy of a non-minimally coupled scalar field, we 
must first obtain the (^-function, i.e. we must perform an integration over 
the fictitious parameter a. Using the formula derived above, the task is to 
evaluate the following integral 

i(s] Jo J 



T{s) Jo 
1 

fis 



(55) 

where we have used the fact that R only depends on the time t and not on any 
of the other coordinates together with the orthonormality of ipxi'x) to carry 
out the integrals over x and the angles. Unfortunately we have not been able 
to perform this cr-integration explicitly, so instead we have to make do with 
an approximation. The expression appearing in the exponential has an easy 
interpretation: the first term is the fiat space-time contribution, the second is 
the non-minimal coupling term while the remaining two terms are first loop 
quantum corrections to the gravitational background. It thus makes sense 
to only expand these latter terms to the first order, whereby equation (^) 
becomes 

A 

-]^i^ies{s + l)(s + 2){iR + A)-^-^) dt (56) 

Whence the effective action becomes 

non— minimal C (0) (^'^) 
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= -E/ /^a^^V^^o^^*' [-HCR + x) + lmR + xy 



(59) 
(60) 



where the curvature scalar and its derivatives are found to be 

R = 6a~^{ad + d^ + K) 
it = -2a~^dR + 6a~^{3ad + aa^^^) 
R = (6a"V - 2a"^a)i?- 24a"^d(3da + aa(^)) + 

6a-'(3d2 + W^)+aa(^)) (61) 

here a^^^ and a^"^^ denotes the third and fourth derivative of a with respect 
to time. With this the pressure, which once more is to enter the Einstein 
equations as (233), becomes 

1 



E 



_1 3/2 dhx 
, 2"' da 



-I 



da 



dt'\ (-ln(^i? + A) 



+\iRm+\)-'-\eR\iR+x)-' 



^2eR^,R^xr'^) 



-2 /J Vf^dt' 



(62) 



The differentiations with respect to the scale factor appearing in this expres- 
sion are 



dR 
da 
dR 

da 
dR 
da 



-2a-^R+6a-^d 



6a dR — 12a dd 



(63) 
(64) 



;i2a~^d^ - 2a~^d)i? + {6a'^d^ - 2a~^d)(-2a~^i? + 6a~^d) 

+72a~^d(3da + aa^^^) - 12a-^(3d^ + 4da(^) + aa^^^) 
-24a-3da(3) + Ga'^a^"^^ (65) 
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and 



dhx 
da 



-6da ^ -\ — a "^a — 2a 
2 



All this is to be inserted into the Einstein equations 

a^^(2aa + a? + K) = p 



(66) 



(67) 



The reader will no doubt appreciate that these equations are not presented 
in their discretised form. 

Again, as with the minimal coupled case, this result differs from the equations 
of motion found by Suen and Anderson []14| due to the use of the WKB- 



approximation, which speeds up convergence. But this time, there is also a 
discrepancy due to the presence of an e~^^" term in the heat kernel, which, 
as demonstrated elsewhere, should be present . 



4.2 Performing the Simulations 

We have carried out the same discretisation procedure as for the minimally 
coupled case and have simulated the evolution of the various Friedman- 
Robertson- Walker geometries for the case of conformal coupling ^ = \- The 
results are shown in figures 2a-c. 

We see that this conformal coupling to the gravitational background makes 
the tendencies already inherent in the minimally coupled case much more 
significant. For K = +1, the closed universe, the tendency towards collapse 
is now even more pronounced, in fact all the chosen values gave rise to a 
collapse, whereas the expansion in the open and fiat universes is seen to be 
faster. 

Again, for large times, we get power-law behaviour (assuming no surprises 
for really large times), giving powers ~ 0.6 for K = Q and 1.5 for K = — 1. 
Note that this means that for a vacuous FRW-geometry with K = —1 and 
with a conformal scalar field as its only inhabitant, the age of the universe is 
no longer bounded by its Hubble time. Also note that for K = our result 
deviates only little from that of Starobinsky (who finds a power of 2/3) 
Let us conclude this section by comparing the evolution for the minimally 
and the conformally coupled scalar field cases, see figure 3, to gain a little 
intuition about the influence of the coupling to the background. 
It is seen that the coupling to the background tends to speed up expan- 
sion/collapse and also that the onset of expansion/collapse occurs earlier. 



19 



Also some some initial conditions that for = (seemingly) give rise to ex- 
pansion in the minimally coupled case gives rise to collapse in the conformal 
case and vice versa some initial conditions that ior K — —1 give rise to col- 
lapse in the minimally coupled case gives rise to expansion in the conformal 
case. 

5 Relating the Casimir Effect of (Spin 1) 
Abelian Gauge Bosons to That of a Scalar 
Field 

In order to be able to actually carry out the computations for higher spins as 
well (and thereby eventually being able to handle realistic field theories such 
as the standard models or various GUTs), we want to derive some relation- 
ships between the heat kernels, and thus the ^-functions, for vector bosons 
and Dirac fermions on the one hand to that of the scalar field case on the 
other. As the case of gauge bosons carries some resemblance to that of the 
non-minimally coupled scalar field we first consider that case. 
We can derive an expression for the effective action for a spin-1 field, even 
in the non- Abelian case, by a simple extension of the technique used for the 
non-minimally coupled scalar field, but we need a way to handle the self- 
interaction (the non- linear, non- Abelian terms in the field strength tensor), 
and we will introduce a way of determining meanfields that allows us to 
handle such non-linearities (conseqeuntly, we could then also handle, say, 0^ 
theory in this manner). 

The self-interaction is treated by a mean field approximation - the mean 
field itself being expressible in terms of heat kernels. Due to the implicit 
complexity of the solution we have not been able to carry out the simula- 
tions for vector fields, however. But it is important to emphasise that this 
trick of using a mean field coupled to the way we compute heat kernels for 
non-minimally coupled scalar fields, actually allows us to find an expression 
for the heat kernel of a Yang-Mills theory, and consequently also the effective 
action of non- Abelian theories. 

In order to keep the notation as simple as possible we have left out the 
ghost contribution coming from the over-counting of degrees of freedom in 
the functional integral. In Lorentz gauge, the ghost contribution is simply 
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—2 times the effective action for a minimally coupled scalar-field. It is thus 
straightforward to include this correction at the end. 

When considering a Yang-Mills field in a curved space-time then, in order 
to obtain the field strength tensor the naive guess is to replace the deriva- 
tives of the Minkowski space field strength tensor with covariant derivatives 
which is not correct as this leads to a non-gauge covariant expression (even 
though, accidentally, it gives the right answer in the case of abelian fields and 
no back-ground torsion). Instead proceed by considering the full theory of 
Dirac fermions interacting minimally with the gauge fields as well as with the 
gravitational field. In order to preserve local gauge and Lorentz covariance 
construct a covariant derivative of the form ^ 

= e^^{d, + ^<(x)X,, + igAl{x)Ta) (68) 

where is the vierbein (local base vectors of a freely falling observer), 
io^^{x) is the spin connection being the gravitational analogue of the gauge 
field v4^(x) and X^q the corresponding Lorentz group (5*0(3,1)) generators 
analogous to the gauge group generators T^. Greek indices refer to curvilinear 
coordinates while Latin indices refer to local Lorentz coordinates (of a freely 
falling observer) and are also used for gauge group indices (it should be clear 
from context: in general we will use small Latin letters from the beginning 
of the alphabet to denote gauge indices, whilst reserving letters from the last 
half of the alphabet for use as Lorentz indices). 

As in fiat space field theory the gauge field strength tensor is obtained 
from the commutator of the covariant derivatives 

[D^, D„] = Sl^{x)D, + '-R^^^^{x)X„ + zF^„T, (69) 

(where S^^{x) is the torsion and R^J"^{x) is the Riemann curvature tensor). 
Using the commutation relations 

[df,, <9J = , alternatively [dm, 9„] = ((9^e^ - (9„e^)ej^(g70) 

[Ta,Tt] = tfabcT" (71) 
[^mn? ^pg] '^"QmpXnq ~\~ 'i'1]npXmq il^nqXmp if]mqXnp ('''^) 

{flmn being the fiat metric) the normalisations 

Tr{TaT,) = (73) 
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and the action of the S0{3, 1) generators upon the latin index m 

-^PQ^m ^Vgfn^^q '^^pm(^q C^^) 

one obtains (after a lengthy calculation) the field strength tensor 

= e^^e:{d,At - d.A^ + igflA^Al) (76) 



The generating functional of the full theory of a fermion interacting with 
a non-abelian gauge field is 

Z = j VA^ j /)^/)^e^Saugcflold+Sfc„nion (JJ) 



Referring back to equation (p8| ) we see that the fermion part of the action 
contain reference to the gauge field so that one a priori cannot carry out 
the two integrations independently. One could probably treat the fermion 
action as a source term when doing the gauge field integration and then 
subsequently do the fermion integration. Instead however, we are going to 
make a mean field approximation to in the fermionic integral, as well as 
in the higher order terms of the bosonic integral (see later). Thus we can 
consider the bosonic and the fermionic parts independently (the meanfield 
in the fermionic path integral is just a function of space-time variables). 
The bosonic part of the generating functional then becomes: 

= J VA^e~^ d''xe'^ei{d^A)^n-dpAlTa-gfabcAlAlT-)ele'^^[dsA-T^-dr,Ad^^^^ 

I)A^exp(-^ j d^x [dmA'^ndpAin - d^AlndqA^Ta - g/defd^A'^nAiA^T^ 
~dnA'^TadpA'^Tf, + dnA'^TadqApTfi + g fd^fdnA'^TaApA'^T^ 
-gUcA^^A'^T%AlT, + gU,AlA'^T%AiT, + g' UM^fA^^A'^T^A^AlTf] rf^^'^^) 

I)A^exp(-|- 1 d^x [AlT^d^i&'A^'' - 9'"A"'')Tfc 

+^/rfe/<A^T^(9M™« - d^A^'')Ta 
+ ^-g'fa,JdefAlA'^T'^A^''A^^Tf] ) (79) 
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Using the commutation relation ( [70D and the normahsation (^) this can be 
rewritten as 



d X 



1 



+gUeaAtAl{d^A^'^ - d^A^'^) + -g'],,J,,^Al^A^''A\A-^ 



0) 



where the last part of the first term can be eliminated by applying the Lorentz 
condition, d^A^ = 0. Now make the following mean field approximation to 
the path integral: 



VA^exp{- 



9 



d X 



Anidm^n — dnG'in)e^ dpA^ — A'^dmO'^A'^ 



+Al{gUUd^A^^ - d'^An)At + A-^{-g'UJ,,fA'i^Ai)A'^'' 



) 



=J VAf,e' / d^xAlsl\^-5-s^dpap+s-id„e 



5.1 Solving the Heat Kernel Equation 

The full heat kernel equation then is 



3_ 

4 



,1 



{^STg'fe,cf,^.A;A^')\ = ~d.G:l^\x,x' ,a) (82) 



or, in short hand notation 

J [s^jT^o - sisrdp - f:^\{A)) 

withQ 



G'n(fe)(a^'a;',a) 



-daG'^l^^{x,x\a) 



(83) 



E^P = (a„e™'^ - d'^e^) el (84) 

Now, in order to be able to make use of the WKB solution from the scalar 
case, we must replace the "flat" d'Alembertian Dq = 9p9^, by the proper 

^ The quantity is the structure coefficients of the algebra of "flat" derivatives dm.- 
Thus it is therefore not surprising that it occurs here. 
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d'Alembertian operator on a scalar field, □. From the definition with 
Qnu = Vab^lj^el and e = ^ being the vierbein determinant, it follows that 



□o = □ - -d^{ee';)dP 



(85) 



Inserting this into 



we get 



-d^G-^^ix^x^a) (86) 



where we have also rescaled a —>■ a = g^a/A, in order to get rid of the over 
all factor of ^ in (|83|) and also introduced E defined by 



1 



Sr = Sr + -^d,{ee^ ^S-) (87) 
The first order term is eliminated, as before, by the substitution 

(88) 



where G is a tensor but the factors on the right side of the equation indepen- 
dently seen are not, because the exponent is not. This leads to the following 
equation 



1 



Gn{l)i.x, x\ a) = (x, x', a) 

(89) 



m(a) , 



which is of the form 



+ 



GnS) (x, x',cr) = (x, x', a) 



(90) 



where O = 0{{A{x)),x) is a matrix-valued function of x. The heat kernel 
becomes (once more) a matrix-valued function which we can assume to be 
of the form 



a 



m{a) 
n{b) 



(x,x',a) = Gn(x,x',a)(e"^(^'^'^)):r(l';) 



where Gn denotes the heat-kernel of □ and T is some matrix (T)^^^^"'' = T^(^f^^ . 
Inserting this expression for the heat-kernel into the heat equation we arrive 
at an equation for T^^"'* 



(91) 

nm{a) 



m(a) \ 
k{c) 



■n{b) ) 



■2apin(Go)a^T:;J;')+a 
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m{a) 
n(b) 



d 



m{a) 



da 



(92) 



where a summation over repeated indices is understood. The third term on 
the left hand side vanishes along the diagonal x = x' hj the argument given 
in section 4. We will furthermore write T^^"'' as a Taylor series 

oo 

^5?(^,^) = E-^^5if(^)^^ (93) 



This leads to a recursion relation for the coefficients t^'^^^^q^^ 

u'=0 

or, as dp = e^dfj, and is the inverse of e^: 

^r^'^:^^ + t (5p-^^-^'lTir^)(5^-^^'S?) = -(- + 1)-^^11? (95) 

u'=0 



Due to the boundary condition one has 

d'r^'^ = l^tSr (96) 
and because of this and equation (R^) the next coefficent becomes 



^(1) ^ ^.g^^rn, _ }_StVr>,£r ^'n' + ( (^) ) (97) 



and subsequently, cf equation (pif), the next few coefficients are seen to be 
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and 



r 



(3) m(a) 
nib) 



2 



'n(6) 



1 



2 



1 . 



2 

Id, 



kp 



1 

2 



2-« 

fVaic) 

" ln{c) 

P •'Kb) ri 

^ pmpr:, n j!k{a) 
4 

-\-fl r) P"^Pf)1 fKa) 
i(a) , 

+ terms of even higher order 

5 f"^ic) ap fk{a) 
'pJk{b) ^ Jn{c) 



-\{d,^rYd'fnib) 



(98) 



(99) 
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and so forth. 

For computational purposes one should, as indicated above, keep only the 
terms that seem most important when doing simulations, i.e., the most im- 
portant pure curvature term, the most important pure gauge field terms and 
the most important mixing/gauge field-curvature coupling termsQ. Thus we 
arrive at the following expression for the heat kernel: 

This formula then expresses the heat kernel for a vector boson in a curved 
space-time in terms of that of a scalar field in fiat space-time and a matrix- 
valued function depending on the curvature. 

One could go on to develop a similar formula for the heat kernel for a spin 
s boson then, instead of a matrix-valued function, we would then get a rank 
2s tensor-valued function. 

5.2 Determining the Mean Field 

In the above we substituted mean fields for the Yang-Mills field in the non- 
abelian terms of the action. We now proceed to determine these mean fields. 

^ Note, however, that by discarding terms that are essentially higher and higher order 
derivatives of the curvature one might lose a large contribution when close to singularities. 
Also note the discussion at the end of this subsection. 
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By definition the gauge mean field squared isQ 

where 5" denotes the appropriate action. We want to carry out this functional 
integral. In order to do this, we note that in n-dimensional Euclidean space 
the following holds[| (see appendix for derivation) 



where M is some symmetric matrix. This result is independent of the dimen- 
sionality n and can hence be carried over to the infinite dimensional Hilbert 



space L'^. Thus, in the continuum case where we have M ~ . c^^A^b ( m ; 



have by inference 

^ 2\6A':^ix)6AUx')J "2 



= I ( n^TT^^TTT-x] = b'^'VmnSix, x')D{x, x') (102) 



with D being the Green's function - (the inverse of the matrix M = |^). 
Now, given an operator A, a simple relationship exists between its heat kernel 
and its corresponding Green's function, namely 

D{x,x') = ~ / Ga{x.,x' ,a)da (103) 

JO 



^Wherever the mean value of a any odd power of the gauge fields occur, we substitute 



{A^) for {A). Of course, for {A) ^ 0, this corresponds to spontaneous symmetry breaking 
of gauge symmetry (at least in the meanfield approximation) due to the fact that virtual 
particles interact with gravity while propagating (see e.g. Q). A couple of things have 
to be said about this fact. Firstly, it has been noted that the effective action need not 
be gauge invariant beyond the 1-loop level Secondly, we do indeed prefer to see the 
symmetry breaking as an indication that a more complete theory would have to have 
some 'mixing' of gauge transformations and coordinate ones, i.e., one should only be able 
consider gauge transformations independently of local Lorentz ones as a 'corner' of the 
theory. Actually, at the higher level one might well find that masses occur naturally in the 
theory, as witnessed by the condensate formation seen in the mean field approximation. 
Thirdly, it is most satisfying that the fact that an observer in curved space do see an 
(Casimir) energy density is related to the fact that mean fields do attain finite, non-zero 
values in curved space vacuum. 

^One should note that this gives an interpretation of the corresponding (^-function at 
s = 1, since we get {x^) = i^A/ (1). 
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We can derive this formula by writing the Green's function (i.e., the left 
inverse of the operator, AD{x,x') = S{x,x')) in terms of the eigenfunctions 
of A: 

D{x,x') = Y.Ulix')M^) (104) 



and then using 



A A 



1 



A 7o 



and the spectral representation of the heat kernel to obtain ( |10tj| ). In order 
to obtain the mean field values {A^A\^ we use, to a first approximation, 
that part of the action which only contains the abelian terms (i.e., excluding 
the terms involving the structure constants, fabc^ of the Lie algebra). If 
one desires a better approximation, the procedure of this paragraph may be 
repeated using the full path integral with the first meanfield approximation 
substituted for the gauge fields. This scheme can then be iterated until the 
desired accuracy is obtained. Towards the end of this subsection, we will 
describe the changes being made in the appropriate expressions during this 
iteration. 



The heat kernel equation is as equation (|89D but with = 0. The solution 



is then (still in matrix notation) to the chosen order 

G(x, x', a) = Gu,{x, x', a)e--^^+^^^'-^''^' (105) 
with (using (0)) the coefficients given by (to the order chosen) 

= (107) 
C = {d,AtWAl) (108) 

For the mean field we then get 
1 

{Al^{x)Ai{x')) = lim-<5'^Vn5(x-x') / dacj'-^G{x,x,a) 

z Jo 

^"^Even though the matrices A, B, C do not a priori commute, there is no problem here 
to the chosen order of accuracy. It is conceivable, though, that one encounters difhculties 
when iterating the scheme for determining the meanfield because the gauge field dependent 
terms are not symmetric. 
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1 

lim -5"'^r]mn^{x — x) / daa' 
s^o 2 JO 



s-2 ^,-Acr+};Ba'^-\Ccr^ 



(47r)-2i5«Vn'5(a; 



lim 



2 



3 



da 



327r2 



X 



^'^'r^^n^lx - x') 

lim f ^-("-i)r(s - 1) + ii3^-(^+i)r(s 



i)-ic^-(^+2)r(s + 2)) 

(109) 



where we have used that A is related to the curvature, and hence we can 
make do with only including the first terms in the Taylor series defining 
exp(|i3(7^ — |Ccr'^), just like we did for the non-minimally coupled scalar field 
in an earlier section. It turns out, however, that the resulting integral is 
divergent in our case (as s — we get a term involving r(— 1), since we do 
not have a l/r(s) factor outside). We regularise this singularity simply by 
using the principal value of the meromorphic function r(z) at z = — 1 (see 
[|]). Laurent expanding we get 



Viz) 



—— + (7 - 1) - -{z + 1)(67' - 127 + vr' + 12) + 
z + I 12 



Since the principal part of a meromorphic function at a pole is just the finite 
part, we finally obtain 



{A'Lix)Ai{x)) 



rcg 



^(7- 

('^mp(7 



1) + ^BA'' 



-CA~ 



1) + -l3mk{A ^)p — -Cmk 



{A^ 



(110) 

(111) 



Now, upon iteration of this we must obviously include the mean field in this. 
We then take the above result for (yl^yl^)i.eg and insert it into the action for 
the Yang-Mills field as an extra term. The heat kernel for this new action is 
then found (the mean fields merely corresponds to an extra curvature term) 
by the same methods as above, from which one gets a better approximation 
to the mean fields, and so on. 
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5.3 Using the method on other problems 

A note on the generality of the above calculation of the mean field is in order: 
We removed the flat space d'Alembertian Dq by introducing the curved space 
d'Alembertian □ through equation (^). This was done to mimic the case of 
the non-minimally coupled scalar field and is possible because we know Gn 
from the section of the minimally coupled scalar field. However, for other 
problems one might not know Guq- It should thus be remarked that had 
we continued using Dq instead of □ the functional form of the heat kernel 
equation (pOD would have been unchanged and the derivations proceed as 
before. One can then take advantage of the fact that the heat kernel of the 
'flat' space-time d'Alembertian, Dq , is known to be (se e.g. [jl], ^, 

Gno{x,x ,a) = _ e « (112) 

in four dimensions. Here A{x, x') is the geodesic distance between the points 
X and x', i.e. A = ds, which in Cartesian coordinates would be simply 
|x — The exact form is not needed, only the fact that Gn^ taken along the 
diagonal x = x' is independent of x. One should also note that one can only 
use freely falling coordinates for as well x as x' when these are sufficiently 
close to each other, i.e., both within some sufficiently small neighbourhood. 
But as we would only be interested, ultimately, in the limit x' —>■ x, this 
approximation would be justified. 

6 Heat Kernel Equation and Zeta-Function 
for Spin 1/2 Fermions 

In the fermionic case the heat kernel equation becomes 

d 

)^G'y(x, x', (t) = ——Gy(x,x',a) (113) 

^^Actually, the formula in jsj is for Cartesian coordinates. Here we have brought it on 
a covariant form. This form is the simplest possible form for the heat kernel compatible 
with the equivalence principle. 
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Now assume that we know the eigenfunctions of ^ i.e. ^{jpxja = ^{"ipxja 
and guess that the heat kernel has the following form 

=Y.c-^{M^')UM^))pe^^'^'' (114) 

A 

which inserted into the heat kernel equation gives 

(115) 

provided that [f, c"'^] = and ^(A) = -A. 

To establish the link with the scalar case make exactly the same exercise for 
the operator i.e. assume that the heat kernel for is 

Gy.{x,x',a) =T.^"^M^')ra{Mx))pe'^'^'' (116) 

A 

and insert this into the heat kernel equation to get (using ^{ipxja = ^{'^x)a 
twice) 

A 

d 

= ——Gy2(x,x',a) if one choose h{x) = — A^117) 
If one choose to normalise the spinor eigenfunctions as follow^ 

E / V9d'x^d'x'c''^{M^')):{4j',{x))p6{x - x') = 6{X - A') (118) 



A 

(if is chosen to be the Kronecker delta then this is the ordinary 
normalisation) then the zeta-functions are related as follows (cf equation 

O); 

Cy2(s) = j daa'^^Y j GY2{x,x,a)y/g(fx 

A 

= Cy(2s) (119) 



choosing another normahsation triviaUy alters equation (117). 
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6.1 Heat Kernel Equation and Zeta- Function for Free 
Spin 1/2 Fermions 

Now note that for a free fermion field the covariant derivative is 

f = re^m{d, + l^7{^)Xpq) (120) 

Representing the 5'0(3, 1) generators in terms of the sigma matrices, Upq — 
I [7^,75], one obtains for the derivative squared 

y2 = (n + ^_^i?).l4 (121) 

(where I4 is the four dimensional unit matrix, and = |) establishing the 
link between the scalar and the fermion cases if one remembers to include a 
factor of 4 (one for each spinor component) in the Dirac case (corresponding 
to taking the trace over the unit matrix). Thus we find 

Cy(«) = Cy^(^«)=4C^+«>(|) (122) 

This is our final result for free Dirac spinors. We have now arrived at formulas 
expressing the heat kernels of free higher spin fields in terms of that of the 
scalar field. We notice, however, that in the fermion case we arrive at a 
non-minimally coupled scalar field equation. The above result is the result 
we will use when performing the simulations. The case of interacting fields 
is discussed in the appendices. 

6.2 Performing the Simulations 

Due to the simple relationship between the (^-functions for a non-minimally 
coupled scalar field and the Dirac field, the simulations of evolution 
of Friedman-Robertson- Walker geometries due to vacuum fiuctuations of 
fermionic fields can be carried out quite easily. As the functional integral 
giving the free energy is now over Grassmannian variables, the sign of the 
energy will be different (we will get the determinant raised to the power of 
plus one half and not minus one half as for bosons), due to the abovemen- 
tioned simple relationship we will furthermore get a factor four (one for each 
spinor component) and a factor 1/2 from the argument of the ^-function. 
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These changes, however, are easily performed and we can carry out the sim- 
ulations with very little change to the program for non-minimally coupled 
scalar fields. We ran the simulations with an initial value of the scale factor 
fixed at 5Lpianck- The results are, for = 0, ±1 plotted in figure 4. Note that 
all = +1 universes either collapse or seem to be on the verge of collapse 
and that the evolution for K = and K = —1 tends to converge. Again, 
this latter part of the expansion is polynomial in time, the corresponding 
power is ~ 0.7. Also note that rather contrary to what one might expect 
due to the fact that the pressure has the opposite sign of that of the case 
of a (conformal or scalar) scalar field (cf formulas (|^ and the corresponding 
one for Grassmann fields for which the power is positive), the behaviour is 
not qualitatively very different in the two cases. Actually, the fermionic case 
looks rather similar to the minimally coupled scalar field. 

7 Conclusion 

We have derived rather simple expressions for the effective actions for quan- 
tum fields of spin zero, one-half and one in a Friedman-Robertson- Walker 
background, thereby obtaining the equations of motion governing the grav- 
itational degrees of freedom i.e. the scale factor a{t) of a vacuous FRW 
universe. Or, in other words, the Einstein equation now describes the cou- 
pling of the zero-point energy (or pressure), the so-called Casimir energy, to 
the Einstein tensor, i.e. we consider the evolution of space-time including 
back-reaction of gravity upon itself as mediated through virtual particles of 
(other) quantum fields. This is not the whole picture (direct gravitational 
self-interaction is of course not included) but probably is as far seeing gravi- 
tation as a background field (first quantisation) can bring us. 
Looking at the plots from the simulations we see first of all that, irrespective 
of whether the residing field is minimally or conformally coupled scalar field 
or a Dirac fermion field, the closed universe K = +1 seems to be ruled out, 
as the simulations show a marked tendency to collapse in this case. 
For K = there is a tendency towards expansion according to a power law 
with power ~ 0.6 — 0.7, whatever the residing field. 

For K = —1 however, the nature of the residing field is rather important: 
In all cases there is a marked tendency, at large times, towards expansion 
according to some power law, the power being ~ 0.7 for fermions, ~ 0.8 for 
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a minimally coupled scalar field and ~ 1.5 for a conformally coupled scalar 
field. 

This indicates two things: First the evolution of even a vacuous space-time 
may depend crucially on its (virtual) content, i.e. on the residing fields. Sec- 
ondly, the quantum field content might conspire to give expansion according 
to power laws with powers greater than 1. Thus, the age of a FRW- universe 
need not be less than its Hubblc-timc. Or stated differently; the local Hub- 
ble expansion could be larger than the one observed in the early universe (at 
large redshifts). 

It should be noted that this gravitational back-reaction through quantum 
fields cannot be done away with, it should be included even in the simplest 
of models that attempts to describe the real universe. We failed to make 
simulations for the case where all fields of the standard model resides in the 
FRW-geometry so we abstain from guessing at the evolution of " a real" FRW 
universe. But we do think that the age-problem of the universe (the problem 
of the large Hubble constant) is non-existent. 

A The Zeta- Function for Spin 1/2 Fermions 
Minimally Coupled to a Gauge Field 

For real world purposes one should consider the case where the fermion field 
couples to a gauge boson field. The gauge invariant and space covariant 
derivative is 

Dm = e(^(9^ + |<(x)Xp, + igAl{x)Ta) (123) 
Using the representation (??) and the relations 

^I'Y = 77^^ + 2z(7f« (124) 

{I'^mn.l'} = 2r^^^am^-r8l5i-25la\-i5l5l-r25la^m (126) 
one obtains, for the derivative squared 
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= □ + e/i? + 2^a^"'F;^r„ + ^V^^p^S7;r6 + 6?(^) (m) 

To make this expression manageable note that putting 

Q{A) = (128) 

is an allowed gauge condition. To show this one needs to demonstrate that 
det(^) 7^ which is easily done. Perform an infinitesimal gauge transfor- 
mation SA'^ = DmUj"- (where cj" is the arbitrary function appearing in the 
transformation rule for fermions ip — > exp{igLj"'Ta)ip in Q to get 

^ = 77™"a^L>" + function 

Since Lorcntz gauge is an allowed gauge condition, we know that det dmD"^ ^ 
0, and adding a function cannot change this fact. Hence Q{.A) = is good 
gauge condition. 

Using the mean field approximation described in appendix A, the gauge field 
dependent terms simply becomes a function which is in principle no different 
from the ^jR term to which it, for calculational purposes, can be added. 



B Calculation of {xiXj) 

We want to derive an expression for {xiXj) where components of 

an A^-dimensional vector, and where 

with M some symmetrical N x A'^-matrix. 

First note that the integral vanishes whenever i ^ j. Secondly, that since M 
is symmetrical we can diagonalise it in which case the problem reduces to 
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the = 1 case. 
For iV = 1 we have 




-1 



(129) 



For arbitrary we then have 



(130) 



where no summation over i is to be performed. This is the needed result. 
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Figure 1: Evolution of the scale- factor for: a)i^ = -|-l, b)ii' = 0, c)i^r = — 1 
with a minimally coupled scalar field. 



Figure 2: Evolution of the scale factor for: a) K — +1, h) K — and c) 
K = —1 with a conformally coupled scalar field 



Figure 3: Evolution of scale factor for the cases of the residing field being 
a minimally coupled (solid line) or a conformally coupled (dotted line): a) 
K = +1, b) K - -1. 



Figure 4: Evolution of the scale factor for : a) K — +1, h) K = 0, c) K — —1 
with a fermion field. 
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